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INTRODUCTION 


In this paper, we extend the analysis of the spectral vanishing viscosity 
method for stabilizing spectral approximations of nonlinear conservation 
laws. The spectral vanishing viscosity has been first introduced in [3], 

OO 

where it was shown that L — bounded spectral-Galerkin approximations converge 
2 

strongly in ^i 0 c^ X,t ^ t0 tke exact en tropy solutions of such conservation 
laws • 

The analysis is performed on the 2TT-periodic inviscid Burgers' equation 


U(x,t) 


,u (x,t). 


submitted to the additional entropy condition 

(1.2) < 0, 

which singles out the unique "physically relevant" weak solution of (1.1) • 
Both the spectral-Galerkin and pseudospectral-collocation methods for (1.1), 

(1.2) are treated, and to this end we proceed as follows. 

Denote by S^u(x,t) the spectral-Fourier projection of u(x,t), 

(1.3) S u(x , t) = l u(k,t)e ikx , u(k, t ) = • / u(x,t)e ikx dx, 

|k|<N 2lf 0 

and let I^u(x,t) denote the pseudospectral-Fourier prolection of u(x,t), 
which interpolate u(x,t) at the 2N+1 equidistant collocation points 
x v = vh, h = , v = 0 , • • • , 2N , 
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2N “-ilex 

(1.4) I u(x,t) = l u(k,t)e lkx , u(k,t) = 77 * I u(x ,t)e V . 

N |k| < N v=0 

These two projection operators differ by aliasing error — that is, we have 

(1.5) 

where the aliasing projection A N is given by [2] 

(1.6) A u(x,t) = T [ l u(k + j(2N + 1), t)]e ikx . 

|k|<N j*0 

Throughout this paper, we use 


(1.7) 


N 


S N + a<A K 


as a concise notation for the two kinds of Fourier projections: either with 

a = 0 , corresponding to the spectral projection, or with a - 1 which corre- 
sponds to the pseudospectral interpolation. 

We approximate the Fourier projection of the exact solution P^u(x,t), by 
an N-trigonometric polynomial, u N (x,t) 

(1.8) u (x,t) = y u,(t)e ikx , 

N |kf<N k 

which is determined by the approximate evolution equation 


at U N (x,t) + 8 x (P N 2 U N e 3x (Q N 3x U N (x>t)) * 


(1.9) 
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The expression on the right-hand side of (1.9) represents the spectral vanish- 
ing viscosity term. Here is the spectral viscosity operator which is 
defined as a convolution with a symmetric viscosity kernel Q N (x), 


(U10) % h V x '° ’ V 10 * fe V x>t) ’ V x) ‘ 


Q(k)e 


ikx 


In the spectral case where a = 0, (1.9) amounts to 


3t u N^ x,t ^ + 3x ^ S N 2 u N^ X,t ^ e 3x ^N^ x ^* 3x u N < ' X,t ^’ 


consisting of a nonlinear system of ordinary differential equations for the 

A 

Fourier coefficients, u^(t), which are coupled through the standard spectral 
convolution treatment of the nonlinear term. The interpretation of the scheme 
(1.9) in the pseudospectral case where a - 1 leads us to 


3t U N^ X v ,t: ^ + Ix^N 2 u N^ x,t ^ |x=x e TFx^N^* 3x u i/ x,t ^ |x=x ’ 0 £ v £ 2N > 

1 V 1 V 


and consists in a complete statement of a standard collocation method with a 
pseudospectral treatment of the nonlinear term. 

In both the spectral and pseudospectral cases, the spectral viscosity 
operator can be efficiently implemented in the Fourier rather than the physi- 
cal space, i.e., 


!„(*). !jU N (K,t)l - -e.|^ 


£ 7J <Q nfc u N (x ' t)) 5 € k [ V x) * fe u N (x - t)1 ’ k2 0(k)u,.(t)e lkx 


An essential ingredient of our spectral viscosity operator, Q^, is that it 
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should operate only on the high portion of the spectrum, in order to retain 
the formal spectral accuracy of the method. Hence we make 

ASSUMPTION I: There exists a constant m = m(N) < ^ N, such that 

| k | < m, 

m |k| jC 2m, 

2m < | k | < N. 


Q(k) = 0, 

0 < < 1 » 
Q(k) = 1, 


Then, with Qjj - I - R,,,, we can rewrite (1.9) as 


( 1 . 11 ) 


Tt u » <, ' t) 


+ — (P — 

3x ' N 2 


u N (x,t)) 


H 1(I - 




where the corresponding kernel R m (x), 


( 1 . 12 ) 


R (x) = I R(k)e ikx , 
® |kj<2m 


is a trigonometric polynomial of degree _< 2m, with Fourier coefficients 


(1.13) 


R(k) = 1, | k | < m, 

,0 < R(k) <^1, m £ l k ! .1 2m* 


In order to guarantee the uniform boundedness of our approximation, u N (x,t), 
we shall need to control the size of this kernel; we therefore make 
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ASSDMPTION II: There exists a constant such that 

(1*14) HR (*)U < Const. logm. 

m L 1 (*> - 

We remark that the assumption of a logarithmic upper bound for the size of 
^(x) is plausible, since typical applications involve R(k) which de- 
crease monotonically to zero and (1.14) is automatically fulfilled in such 
cases; consult Appendix A. To obtain, with the help of Assumption II, the 
promised uniform bound on ujj(x,t), necessitates L°-bounded initial data, 
uj,j(x,0). For technical reasons we shall need the slightly stronger 

ASSUMPTION III: There exists a constant such that 

“ u N (x,t=°)H < Y |u (t=0) | < Const . 

N L (x) | k J <N k - 0 

The spectral viscosity term on the right of (1.11) depends on two free 
parameters: the viscosity amplitude e = e(N) and the effective size of 

the inviscid spectrum m = m(N). These two parameters should be chosen to en- 
sure the convergence of the method. In [3] it has been proved that in the ab- 
sence of such viscosity term, e = 0, strong as well as weak convergence to 
the exact entropy solution fails. 

The main result of this paper asserts 

Theorem 1.1: Consider the Fourier approximation (1.11) of either 

spectral or pseudospectral type. Let the spectral viscosity in (1.12) - 


(1.14) be parameterized with (e,m) as follows 
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(1.15) e = e(m) 


a 


m 


m(N) 


Const .N 


3 


0 < 6 < 


m • HR (OH 1 
" L l 00 


Then u^(x,t) converges boundedly a.e. to the unique entropy solution of the 
conservation law (1.1). 


Let us examine for example the viscosity operator = I - S m . Here 

R^x) coincides with Dirichlet kernel D m (x), where [5, Chapter II] 


d (x) = l 
m |kf<m 


ikx 1 sin(m + 2 )x 


2tt 


sm y x 


ID (OH t 9 logm, 

m T O \ 1 

L (x) 7 T 


so that Assumption II is fulfilled and Theorem 1.1 yields 


Corollary 1.2: Consider the Fourier approximation 


(i.i 6 ) 3t u N (x,t) + 9x (P N 2 u N (x,t)) e 9x [(I - S m ) 9x u N (x,t)], 


with 


N 


-23 


(1.17) e = e(N) ~ Const. y Q - ^ » m = m(N) ~ Const. N , 0 < 3 < 


Then u^(x,t) converges boundedly a.e. to the unique entropy solution of the 
conservation law ( 1 . 1 ). 


The spectral portion of this result (a = 0), was derived in [3, Theorem 
4 . 1 ] under the assumption that the numerical solution ujj(x,t) remains uni- 
formly bounded. The extension of Corollary 1.2 includes the pseudospectral 
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approximation (a = 1), and in addition, thanks to the slightly more stringent 

parametrization than that of [3, Theorem 4.1], contains a proof of the previ- 
00 

ously assumed L -bound. 

A 

In the last example the viscosity symbols Q(k) were discontinuous 
at | k | = m. It was suggested in [3] that the use of viscosity operators 
Qjq with smoothly varying symbols would be advantageous for the spectral vis- 
cosity method in (1.9). As our second and final example we consider the 
simplest viscosity operator of this type, namely 

I k | < m, 

m _< | k | _< 2m, 

2m < | k | _< N. 


Q(k) = 0, 

Q(k) = , 


m 


Q(k) = 1, 


This kind of spectral viscosity is intimately related to the Fejer operator 
1 m-1 

F m * — • I if we let ^ m (x) denote the corresponding Fejer kernel [3, 

k=0 

Chapter III] 


K_(x) = I 

| k j <m 


m 


(1 ~ 


[k|_ )e lkx 

m 


. 1 

0 sin -rr mx 2 

= (-r4-) • 

sin x 


IK (•)! , = 1, 

" L*(x) 


then for Q N = I - we have R^x) = 2K2 m (x) - K^Cx). Hence the kernel 

associated with 

2m- 1 

R = 2F 0 - F = - • T S 

m 2m m m l 


m m . ^ k 
k=m 


is L^-uniformly bounded, 


11 k ( • ) J! 

” l‘(x) 


< 2,K 2» ( - ), i 1. , 

L (x) 


UK (•)« < 3, 

■ ihx) “ 
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so that Assumption II is fulfilled and Theorem 1.1 yields 


Corollary 1.3: Consider the Fourier approximation 


2 m _ j 

( 1 * 18 > ft U N^ X,t) + 37 ^ P N 2 " 6 9x ” m \ ^ V 3x u N (x,t ^’ 


(1.19) e = e(N) ~ Const. N m = m(N) ~ Const. N®, 0 < 3 < ^ 


Then u N (x,t) converges boundedly a.e. to the unique entropy solution of the 
conservation law (1.1). 


The paper is organized as follows. In Section 2 we derive a couple of 
basic L^-type a'priori energy estimates. In Section 3, these estimates are 
used in order to study the spectral decay rate of the Fourier coefficients. 
This in turn enables us, in Section 4, to obtain L°° a priori estimate on 
the numerical solution. Finally, based on the a priori estimates prepared in 
Sections 2, 3, and 4, Theorem 1.1 is proved in Section 5 along the lines of 
[3], using compensated compactness arguments. 


2. I/-TYPE A PRIORI ESTIMATES 

We consider the approximate Fourier method (1.9) which we rewrite as 


(2,1) at U N + 3x ( 2 U N ) ” 3x [(I V 2 U N ] + E 3x tQ N 3x U N^ E 1 + II# 
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In order to prove the convergence of this method we need a couple of a priori 
estimates on its solution. To this end, we multiply (2.1) by u N 


3 A 2 
3t ( 2 U N 


a_ A 3, 

9x 


<4 O + (y Uj) 


( 2 . 2 ) 




„ . 1 2. 

P N 2 U N 


+ eu 


N 3x 


[Q N 35 U N> 


= III + IV, 


and integrate over a 2ir period: the integral of the second term on the left 
vanishes by periodicity, and after integration by parts for the second term on 
the right we are left with 


(2 • 3 ) 

THt" u N ( *’t)» 2 2 + E / !x V*» t)Q N !x V x » t)dx = l U N ’5x"^ (I-P N ) T u N ]dx * 

L (x) 0 0 


Using (1.7) and the fact that I - S N is orthogonal to our N-space, we find 
that the right-hand side of (2.3) equals 


l “ntj l(I - p N> •z u N ldx * **• I 


and by the aliasing relation (1.6), this does not exceed 


l u NTj [a -Vr u N ldx 


a 

1 


v 

f AAA 

JLmJ IPU U U < 

|p+q+rf=2N+l P ^ 


1|12 2 

|p+q+r|=2N+l 



In view of |p + q + r| = 2N + 1, at least two of the three indices |p| < N, 
f q | _< N and |r| <_ N are greater in absolute value than y , 


and hence 
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t 0 \ >>* 


<¥ £ £ wi: 

|<|p|<N |<|q|<N 

1/2 


A A 


q K Mu 


q 1 1 ±(2N+1 )-(p+q) 


<4^-[ ]C p2|u p |2 S q2,U q |2] ^ 2 ,U ±(2N+l)-(p+q) 12 ^ 

N — - |<|q!<N |<|p|<N |<|q|<N 


1/2 


fe|p|< N 


2^i F t~' 2- 

Consequently, since Q(p) = Q(q) = 1 for |p|,|q| > | , the expression on 
the right of (2.3) can be upper bounded by 


< 2 '*‘>/\ fc' I - P N ) I ^ - |IQ N h 


Moreover, since 0 K Q(k) <( 1 we have for the second term on the left of 
(2.3) 


2ir a g 

(2.4b) ef V X ’ t)Q N 3^ U N (x,t)dx = 


A A 


Y k Q(k)|u k (t)| > elQ N u n (* .t) 11 2 ' 

c | <N h \ x ) 


Inserting this together with (2.4a) into (2.3) we end up with 


2 

(2 - 5) 7 3T% ( -' t), ^ (x) + [* ■?7T ,u N ( ’ ,t>, l 2 (x) 1 ' ,Q N^ U N ( ‘’ e), L 2 (x r 


Thus, as long as 


( 2 . 6 ) 


e - 


72 " u N (,>t)l T 2 

Li 


(x) 



y 


we obtain 
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(2.7) 


d 

It 


L (x) 


+ ellQ N U N ( * >t)l|2 2. . 1 °* 
L (x) 


In particular, (2.7) implies that for 

O 

our first L^-type a'priori estimate 


u M (x,t) = l u. (t)e 

N |k|<N k 


ikx 


we have 


(2.8) lln(.,t)ll 2 _ = l | u. ( t ) | 2 < E 2 , E = lu M (.,t=0)H <Const n . 

N l/(x) |k|<N k U IT (x) U 

Hence (2.6), (2.7) and consequently (2.8) prevail for all time provided (2.6) 
is valid at t * 0, i.e., we require that in the pseudospectral case where a 
= 1 we shall have 


(2.9) 


e(N) > 8E q »N 1/2 


indeed. Assumption II tells us that this requirement is fulfilled, at least 
for sufficiently large N, for 


( 2 . 10 ) 


N~ 2 e -1/2 1 

£ > ConSt ’ TogN > 8 ’V N ’ 20 <1 


Furthermore, temporal integration of (2.7) then gives us the second a-priori 
estimate 


(2.U, <s>t) - « I < * 


Jq < Constp. 
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3. THE DECAY RATE OF THE FOURIER COEFFICIENTS 

Our Fourier approximation (2.1) 


(3.1) 


at u n 


a_ 

3x 


,1 2 n - 3 r /•t 
^2 U N ax 


P N ) 2 U N ] 


+ e 


a 

ax 


3x ’ 


3 12 

consists of two kinds of errors. The first term, I = y^ [(I - P N ) y u N ], 

represents the discretization error, which includes spectral truncation 
9 12 

errors y^- [(I - S^) y u^] as well as additional aliasing errors 
3 12 

-a»y^- [A^ y u^j] in the pseudospectral case. In this section, we borrow from 
Kreiss [1], in order to show that due to the second error term of spectral 
vanishing viscosity, II = e y^ [Q^ y^ u^] , there is spectral decay of the 
Fourier coefficients | u^, ( t ) | , |k| > y N, and therefore, that the discretiza- 
tion error is spectrally small. 

We begin by taking the I - S 2 ^ projection of (1.9): for k > m we have 
by Assumption I, (I - S2k)QN = I ~ S2k> and hence 


0.2) S 2k > u N ] + ( I s 2k ) fe t P N 7 U n] " e fe[ (I ~ S 2k ) T3E U J ’ m < k < N. 


Multiplying by (I - S2k)u N and integrating by parts over a 2it-period, we 
find that 


1 11(1 " S 2k )u N ( *’ t)l ' 2 2 r . 

L (x) 


2ir 


7 l ^ I_S 2kVx “n'^ 1 S 2k^ P N U N dx 6,1 ^ S 2kVx U N^ * ’ ^ 1 2. ** 
0 L (x) 


(3.3) 
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The first integral on the right does not exceed 


7 [ ^ S 2k^ 7x U N* ^ S 2k^ P N U N dx - 


(3.4) 


^I“ (I " S 2k>feV-’ t) ' T 2, / ,(I - S 2k )P N U N ( ’’ t)B 2 . ' 

L (x) L (x) 


In order to estimate the second term of the last product, we will make use of 
the following lemma whose proof is postponed to the end of this section. 


Lemma 3.1: Let f^ = f^(x) and = g^(x) be a couple of N- 

trigonometric polynomials. Then, for any 0 < 2k < N we have 


11 ^ S 2k^ P N^ f N g N^' L 2^ x ^ - 

(3.5) 

[,f " , L 2 (x) " <I - S k4 g N' l 2 (l) + , *N , L 2 (][) -'< I - S k )4E f »' L 2( 1!) I 


Lemma 3*1 with f_ T (*) = g XT (*) = u__(*,t) implies 

N N N 


(3.6) I (I-S 2k )P [) „2(., t )i 2 , , <^.u N (.,t). •'«-S k )| ; u N (.,t) 1 

L (x) L (x) 


L (x) /k 


Equipped with (3.4), (3.6), and (2.8) we return to (3.3) to find that 


7 IF s + 

L (x) L (x) 


+ ^ V ,(I - S 2k> lx ,•'« - V lx V-^',2, . 

» k L (x) L (x) 


(3.7) 
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which brings us to 


Theoren 3*25 For any Integer s _> 0 there exists a constant C g = 
Const(s,Eo)» such that for sufficiently large N, N > 2 s *4m, we have 


Proof : Let E k (t) abbreviates the quantity 


(3.9) 


\ M 5 1(1 - S k )u N (, ’ t), I 2, . 

L (x) 


In view of the inverse inequalities 


2kE (t) < .(I - S ) ^ UjjC* »t) ' . < NE 2k CO, 

L (x) 


it follows from (3.7) that E w (t) satisfy 


j 7 2E n -N 2 

<3 ' I0) « E 2k (t) < - 4Ek E 2k (t) + \ M ’ 


m < k < N. 


Temporal integration yields that for any 0 < tg ^ c we have 


2E r »N 2 t . .2, . -4ek 2 (t-t ) 

(3.11) E 2k (t) < — / e" 4ek (t T) E k (x)dx + e 0 ‘E^t^, 

^ T=t 0 


and therefore 


E„.N 2 -4ek 2 (t-t n ) 

(3.12) E 2k^ fc ^ Z ~ 7 * max E k^ T ^ + 6 * E 2k^ t 0^ 

2e/k*k tQ<x<t 
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The a'priori estimate (2.8) implies that 


max E (x) < max E (x) < E , 
0<T <t ^ 0<T <t ° ” ° 


and in view of (3.12) we have 


(3. 13) 2k 


E n* N —L 

E 2k (t) i 7"/=''2 * E 0 + e 6 * E 0> k > m * 

2e/k»k 


If we choose t Q - in (3.12) we find that 


(3.14) E. (t) < 


V N 


-2ek t „ ,t, 


2k' w “ 2e/k» k 2 * t^ Ek(T) + 6 ~~ fc ’ E 2k ( 7 ) » k > 


and following Kreiss [1], we can use this to improve our estimate (3.13), 
Namely, for k > 2m we can use (3.13)^ to upper bound the terms max e (t) 

K 

t T 

an< * on t ^ ie (3.14), and obtain the improved bound 


8E f,*N 2 2 8E n .N 2 ,2 

E 2k (t) < • E 0 + C 1 + -r-^) e e k > 2m. 

s/k»k 


:/k*k' 


Now we can repeat this process, and by induction we obtain that for k > 2 s »m 
we have 


8 s E rt *N 2 s+1 


w “n " 8 S E-»N 2 s .-s+1 , 2 

(3 - 15 > - Eo + (l + — ***..„. 

e/Tc*k e/Ic*k u 


Verification of the induction step is left to the Appendix. Finally, (3.15) 
implies that for sufficiently large k = ^ > 2 s »m, we have 
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1(1 - S N )u N (.,t)ll 2 = E N (t) < 

T L (x) -r 


32* 8 S *E n s+1 32»8 S »E n s ,-(s+l) 

< ( it-—) -e„ ♦ (i + «' 4 ‘ £N ‘ 

e/N e/N 


V 


and 3.8 follows. 

Our parametrization in ( 1 . 14) , ( 1 . 15) implies that for sufficiently 
large N we have 


4 ^ S+ ^*eN^t > Const. 


2 ~ 2 0 

N _ . t > N 3/2 .t, 


4 S+ ^ *logN 


0 < 2B < 4 , 


as well as 


< Const. — * < — , 0<Y<i"28, 

e/ N “ d -26) - N y 7 

N 

and Theorem 3. 2 tells us that 

Corollary 3*3* For any Integer s > 0 there exists a constant C g 
such that 


3/2 

(3.16) 1(1 - S N )u (»,t)l 2 < C s *(N _S + e“ N * fc ). 

Y l/(x) S 

Corollary 3.3 indicates the spectral decay of the Fourier coefficients 
|u(k)| with wavenumbers |k| > which in turn implies a similar decay 

for the discretization error, I, on the right of (3.1). For the latter we 
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have 


(3.17) l(H II 4»J('.t)l , * la-s^ujc.,')' 2 , + a 2 .A N ^u 2 (.,t). 2 . 

L (x) L (x) L (x) 


The Fourier coefficients of the two expressions on the right are given 
respectively by 


((I - VT U N (, »t»k“i 2 u (tOu (t) , I k | > N, 

p+q-k=0 y H 


1 u 2 (* ,t)) k =1 2J u^(t)u^(t). 


^ 2 U N 


| p+q-k =2N+1 


P q 


| k 1 < N. 


In both cases, either |p| > ^ or |q|>^; hence each one of these 
coefficients can be upper bounded in a standard fashion to yield 


1 2 


2 ... 1 2 . 


C 1 s N )-j U N (* 2 . v +a IIA N7 V' 1 ^ 1 2, .1 
L (x) L (x) 


y (l+a 2 ).[ y |u(t)| 2 ].[ I I u (t) I ] < 4E 2 .N«ll(I-S )u (*,t)ll 
N<|k|<2N |q|<N q |p|> N “ Y L 2 (x) 

-—2 

and by (3.17) this is the same as 


(3.18) ll(I-P N )| u 2 (*,t)ll 2 < 2E q ./F *1(1 - S N )u N (*,t)» 2 . 

L (x) j L (x) 

Corollary 3.3 together with 3.18 show that due to the presence of the spectral 
viscosity terra II on the right of (2.1), the discretization error I decays to 
zero at a spectral rate independently whether the underlying solution is 


smooth or not. We state this as 
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Corollary 3.4: For any integer s 0 there exists a constant C g 

such that for sufficiently large N we have 

3/2 

(3.19) II = f(I-P N 4 uj(»,t)]l 2 < C g ./N .(N~ S + e" N * t ), s > 0. 

L (x) 

We close this section with the promised 

Proof (of Lemma 3.1) : Starting with the identity 

5 f N (I -V% + (H k>VW Vn-Vh 

and subtracting from this (I - S 2i c ^ P N^ S k f N* S k g N^ = (I “ S 2k^ S k f N* S k g N^ = °» 
we can write 

(I - S 2k )P » (f N S N > 5 

The quantity inside the right brackets is a trigonometric polynomial of degree 
< 2N and hence, by Parseval relation, its L^(x) norm dominates the L^(x) 
norm of its projection, i.e., 

y(i-s 9 . )p (f g )i 9 < ip t...]i ? < 

2k N N N l 2 (x) - N l 2 (x) _ 

(3.20) 

< 'V^VV (I -VV s k%\2, • 

L (x) 


The norm on the right of (3.20) is upper bounded by 
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(3.21) 


" f N ( ' I " ::> k^ g N +(I_S k^ f N* S k g N 11 2. S 

L (x) 


1 f N " 2 . * ,(I - S k )g N" », + " g N " 2 , / ll(I -W», / 

L (x) L (x) L (x) L (x) 


Finally, for hjj equals either fj^ or g^ we have 


(3.22) 


V (I-S, )h„l < ) h < 

k ” L*°(x) - |pt>k P 1 - 


<[. I \] 1/2 -[. I P 2 |J | 2 ] 1/2 <i; KI-S 4 V 
| P T>k P |pj>k P /k k9x N l 2 ( 


and (3.5) follows from (3.20), (3.21), and (3.22). 


4. L A PRIORI ESTIMATE 

The classical energy method can be used to show that the solution of 
(2.1) remains uniformly bounded during a small finite time. The method 

reflects the fact that for sufficiently smooth initial data, say with 

g2 9 

— s- u (x,t=0) which are L 2 -bounded, the process of shock formation takes a 

9x . a 

finite time, during which u^(x,t) remains uniformly bounded and a 

couple of Sobolev norms could be a'priori estimated during that time. 

For a brief initial time intervals, we can do better with regard to the 
smoothness of the initial data, as told by 


Lemma 4.1: Consider the Fourier approximation (1.11) - (1.14) with 

initial data u«(x,t=0) such that Assumption III holds, i.e.. 
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(4.1) 


I |u k (t 

c|<N 


=0) | £ Constp. 


Then for 


t < — we have 
~ N 


(4.2) II 


U » ( ’’ t), L-( X )- 2 '|4 |Uk(C " 0)l ' — 8Const 0 »N 


Proof: The Fourier transform of (1.9) reads 


(4.3) 

u k (t) + ik t + ^_ k u p (t)u q (t) + a * 


E 


u ( t) u (t)] = -ek Q(k)u (t). 
|pfq-k|=2N+l P q k 


Multiply the real (imaginary) part of this by sgn(Reu k (t)) (respectively, 
sgn(Imu. (t))) and sum over all k's: since the right-hand side is negative we 

obtain after summing both parts 


(4.4) 


|r . I < (1 + | a | ) • 2N • T |u (t)|* l |u fc (t)| < 

dt | kj <N k |p|<N p k K p 


< 4N *(. l u k (t)| ) ♦ 


I k l<N 


Integration in time yields 


i u..( • , t ) i <_ y |u,(t)| < s • y | u, ( t=o ) i , 

N L (x) | k] <N k l-4Nf I | u ( t=0) | |k|<N k 

| k I <N * 


and (4.2) follows. 
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To obtain L bound for later time, we shall carefully iterate on 

the L,P(x) norms of u^(x,t). To this end, we multiply (2.1) by pu^ ^ 

and integrate over the 2iT-period, obtaining 


<L, 

dt 




(4.5) 


L F (x) 


■ rf “ST 1 fe [(I ‘ P N ) 7 u N ,dx + p 'l u n"‘ £ k [q n u N |dx - 


2tt 


12,,. r P-1 3 r. 3 


By Corollary 3.4, the discretization error is negligibly small: using (3.19) 

1 2 

and the fact that (I - P ) ^ is a trigonometric polynomial of degree 

< 2N, we have for any s > 0 


(4.6) < 

1/ (x) 

3/2 

™ -'fe 1<I - P N ) 7'4 ( -- t) " I 2, < C s N 2 .(N- s + e- N ' t ). 

L (x) 


Therefore, by Holder inequality, the first integral on the right of (4.5) does 
not exceed 


2n 


P 'l “ST 1 fe I(I - P N , 7 u N ,dX < ! ” lu r 1< ‘ ,t), L q (x )' l |xt <I — P N ) 7 U N ( '’ C) l l i p (!i) £ 
(4.7) 


< p*llu XT (»,t)ll P " 1 »C *(N 2 S +N 2 e 
N L P (x) S 


2-s.„2_-N 3/2 .t I + I. 

p q 


The second integral on the right of (4.5), with = I - R,,,, equals 


2ir 


p */ U N le ‘51E C(I " R m ) lx J N ]dx " 


2 tt _ 9 3u„ 2 tt , -2 

-ep(p-l)-/ u p - (^-) dx - ep-J uj" — j(R m u N )dx. 
0 0 3x 
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The first terra on the right hand side is negative for any even integer 
p > 2; for the second term we use Holder inequality as before, obtaining 

2tt 2 2 

(4.8a) ep.J ' u^ 1 <* m ydx < 

Now, since R u = R (x)*ir T (x,t) is a trigonometric polynomial of degree 
9 m N m N 

< 2m, (consult (1.12)), we can estimate the L?(x) norm of its derivatives 
as follows [5, Chapter X] 


(4.8b) 


II — r-(Ru) 
3x 2 m N 


< (2m) 2 * HR ( • ) ^*u,, ( • , t ) H < 

lP(x)“ m N L P (x)- 

< 4m 2 IIR (*)ll . - II ul t ( * , t ) « „ 

m l/(x) ^ L P (x) 


Using (4.8a) and (4.8b) we conclude that the second integral on the right of 
(4.5) is upper bounded by 


(4.9) p./ 2 ' [Q„ u N )dx < 4p.em 2 .|R ni (.)ll Ll( ^-lu N (..t)lP p( ^. 


We recall that according to our parametrization (1.15), enr*IIR (*)ll , < a. 

” L L (x) 

Hence, equipped with (4.7) and (4.9) we return to (4.5) to find that 




3/2 

pllu N (*,t)ll p ~* [4a»lu N (»,t)l + C g *(N 2-S + N 2 e" N * t )], 


or, after division by the common factor of p II u^( • , t ) II p * , 

S L P (x) 
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d_ 

dt 




L P (x) 


< 4a 


«u N (. 


,t ) II 


+ C 


l/(x) 


, 2-s 2 -N 3/2 .t. 

(N +N e ). 


Finally, we integrate in time obtaining by Gronwall's inequality that for 
any 0 < t Q < t. 


(4.10) 


4a( t-t Q ) 

e • 


, v-^o ) " L p (x ; 


c 

s 


(N 2 S .(t-t Q )+ /N 




Letting p even tends to infinity, then (4.10) with 
gives us 


t 0 t 0 (N) 8Constp»N 


(4.11) lu (• ,t) I < e 4at .[llu. T (*,t_)ll + C (N 2 S *t + /IT e _Const * /f )l 

N T °°, * — L N U go s 7 J 

L (x) L (x) 


and together with Lemma 4.1 we conclude with the desired l“ bound, namely, 


Theorem 4.2: Consider the Fourier approximation (1.11 )— ( 1 . 14 ) . Then for 

any s > 0 there exist constants a > 0 and C such that 

— S 


(4.12) 

%(-,t)» „ < e 4at .[4. J I^Ct-0)! + C *(N _S .t + /TT e-C 0118 *-^)]. 

L (x) |kj<N ° 

Remarks : 1. We observe that the exponential time growth in (4.12) does 

o 

not exceed 4a where a ~ em • HR (• ) II . < Const. 

m T 1 / \ ~~ 

L (x) 
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2. The a'priori LP(x) estimate derived in (4.10) is valid for arbitrary L 2 - 

00 

bound initial data. We note, however, that the resulting L bound in such 
case is not uniform with respect to the initial time tg« That is, with arbi- 
trary L^-bounded initial data, the solution t1u^(*,t)!l ^ may still grow 

L (x) 

like 0(/N). The point made in Lemma 4.1 was that with slightly strengthened 
regularity assumption on the initial data 

A 

I | u ( t=0 ) | < Const- , 

I k I < N 

this growth is bounded for a brief time interval of length ~ ~ , after which 

00 

the spectral viscosity becomes effective and guarantees the L bound later 
on. 


5. CONVERGENCE TO THE ENTROPY SOLUTION 

We follow [3], using compensated compactness arguments to conclude that 
u N (x,t) converges to the entropy solution of (1.1), (1.2). 


Proof : (of Theorem 1.1). Consider the four terms on the right-hand side 
of (2.1) and (2.2). Using (3.18) together with (2.11) along the lines of [3, 
Lemma 3.1], we find that the term I satisfies 


(5.1) 


11 1 5J [(I - p N4'V < -> t)1 ' -1 , 

H loc (x ’ t) 

< ll(I-P„)-i4(.,t)» „ < 4E n .J n . . 

NT ^ L? (x,t) - 0 0 /eWFW 

loc 


0 . 


Also, by the a'priori estimate (2.11) we have 
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(5.2) 


111 5 e h [0 N h U N ] " -1 . 

H l 0C (x > t) 


- ellQ N U n" t 2 


< /eTn7 • j 


L loc (x,t> 


0 


0. 


N-h» 


Thus, in view of (5.1) and (5.2), the terms I and II on the right of (2.1) be- 
long to the compact of H^ c (x,t). 

Next we note that since Qn + Rm = I we have 

/F ll-l— u.J 0 < /F IIQ M 4- u M ll „ + /F III- (R u„)l , 

L loc <x > t) L loc <x>t> L loc ( ,t> 

The first term on the right is bounded by Jq; the second one — being the 
derivative of a trigonometric polynomial of degree 2m, does not exceed 
/e" • 2m •IIilJ 0 < Const. Consequently, 


(5.3) 


/ 'bxV.2 , . i Const ' 




Equipped with (5.3) we now turn to consider the right-hand side of (2.2). For 
the third term in (2.2), 


(5.4a) 


- Vt U N ] 


3 r /t t, x i 2. a 
Tx [u N (I " V 2 U N ] " 9F U N ,(I 


V 7 “n 5 m i + m 2 - 


we have by Theorem 4.2 and the estimate (5.1) 


-26- 


,m i 5 h [ > a " P n4 (x t) - 

(5.4b) l0C ’ 

< -vi-v^v Cx>t) - " u » ( ’’ t>, i."( X r ,<I 

loc 


P)iull 2 *• 

N 2 N L. 2 (x,t) N-*>° 

loc 


o. 


and together with (5.3) we also have 


,m 2 5 k U N • (I ' Vi ’4' 1 . t) ^ 

L loc U ’ t; 

(5.4c) 

(x , t) -^ (I - p ^^ oc ( X , t > 

loc J- oc 


< Const. 4E J n * ** 

“ u e(N)*/N N-*° 


Finally, for the fourth term in (2.2), 


(5.5a) 


IV E eU N fe (Q N hi V 


e h tu N Q N h U N ] 


e lx U N* Q N I x U N E IV 1 + IV 2’ 


we 


have by (2.8), (2.11) and the uniform bound in Theorem 4.2, 


(5.5b) I IT, e £ 7x [U N^ 7x U n" h -1 (* t) - "N 

loc 9 


< II u XT ( - , t ) II */e(M) • J n — •+ o, 
L (x) N-hx. 


while IV 2 = -e [Q N U N 1 2 + £ (I U N* Q N '5x U 1 


satisfies 


■I V 2> 1 , <‘ ,Q n!;V i J 2 , ,v 

L loc (x ' t> L loc <x ’ 5 


+ ^ IIR m 3x ^"2 


^® ,Q N 55T V,2 

J loc' 


Lj oc (*,t> " " W_(x.t) 


(5.5c) 
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+ /e 


•2iii*Eq»Jq _< Const. 


Therefore, by Murat^s lemma, [4], the inequalities (5.4), (5.5) imply that the 

terms III and IV are also in the compact of H^ c (x,t). In summary, we have 

shown that the right hand sides of ( 2 . 1 ), ( 2 . 2 ) lie in the compact of 

H^ c (x,t), and, according to Theorem 4.2, that llu^(*,t)ll ^ is bounded 

OC * L (x) 

(in fact, II u II with p > 6 will do for our purpose). Hence we can 

L P (x,t) 

apply the div-curl lemma [4] to the left-hand sides of (2.1), (2.2) and obtain 

o 

that (a subsequence of) u N (x,t) converges strongly in L loc ( x > t ) t0 a wea ^ 
limit solution *u(x,t). 

Moreover, we claim that this limit is the entropy solution of (1.1). To 
verify this claim we show that the right-hand side of (2.2), III + IV, tends 
weakly to a negative measure. Indeed, by (5.4) and (5.5b) the terms III and 
IV i tend weakly to zero, and hence it is left to show that the term IV 2 > 


IV 


2 


= -e 


V Q N U N 


= -e(Q* 


R ) 

m 


3 3 

V q n JZ V 


N 


R = 
m 


I, 


tends weakly to a negative measure. To this end we proceed as in [3, Section 
4] and rewrite IV 2 in the form 


(5.6) IV 2 - e[Q N ^ u/ e (W’W + e 

3x 


^m'V’VV 


Denote the three terms on the right of (5.6) by IV 21 , ^ v 22 an< ^ IV 23 > 

respectively. By (2.11), IV 2 ^ tends weakly to a negative measure 


N+°° 


(5.7a) 
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If we integrate the second term, IV 2 2> against any Cq test function 
4>(x,t), we find 


/ J +-IV dxdt < * / / l|| • |j <VN ,-q N U Nl d,tdt i 

X t X t 


< e, H" - 

L loc 


, . tl - ,Q rt " L 2 (xt) > 

(X, t ) L ioc (x,t) L loc (x ’ t} 


and since is a trigonometric polynomial of degree <, 2m, this is less 


than 


(5.7b) / / ij;.IV„dxdt < e*ll|^ll .2m* Hu IT 

it 22 " L l0 c (x>t) L loc (x ’ c) 


K. Const. — 7Trt * 0* 

N+00 


Finally, for the third term 


2 g 2 

IV 23 " e — f ^ R m u N^* Q N^ U N -U ^ + 6 TT = IV 231 + IV 232 

3x 9x 


we have 


(5.7c) «IV 231 s e — y ( R m u N ),Q N ( V u) " L 1 ( t) - 

loc * 


3x 


_< e*4tn # l! R m u N ll 2 -llt^-uH 2 ^ < Const. II - u,, 2 


+ 0 , 


\Jl (x,t) L (x,t) 

loc loc 


Hoc '*’ 0 


and since e R m u N *Q N u tends weakly to zero, so does the term ’^232' 


(5.7d) 


wlim [IV 232 = e — y (R m u N ).Q N TT] 0. 
N+*° 3x 
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From (5.7a) - (5.7d) we conclude that the term IV 2 In (5.6) — and therefore 
that the right-hand side of (2.2), tends weakly to a negative measure. Thus, 

by taking the weak limit of (2.2) we recover (1.2) for our limit solution 

— 2 — 

u(x,t). Consequently, the strong L 1qc limit of u N (x,t) = u(x,t) is 

the unique entropy solution of (1.1) as asserted. 
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APPENDIX 


A. THE L 1 -LOGARITHMIC BOUND OF MONOTONE VISCOSITY KERNELS 


We consider symmetric viscosity kernels of the form 


V x) ■ 


Q(k)e 


ikx 


kT <2m 


2ro<fk| <N 


ikx 


with monotonically increasing Fourier coefficients. Then, the kernels which 
correspond to R m = I - Q N , are symmetric polynomials of degree 2m 


(a . 1 ) 


R (x) = 2 
m 


|k| <2m 


R(k)coskx 


whose Fourier coefficients are monotonically decreasing, compare (1.13), 


(a. 2) 


1 > R(k) + > 0. 


Such kernels satisfy Assumption II above, as told by 


Lemma A. 1: There exists a constant such that 


(a. 3) 


HR (•)! , < Const. logm. 

m T I / v 
L (x) 


Proof ; The result follows if we can show that R^x) is majorized by 


Const .m and Const, , for then we have 
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(a.4) 


(•)H . < / Const .m»dx + / Const, -p-r dx _< 

" i< w <, 1x1 


x<- 


_< — Const .m + 2 ‘Const. log IxJl^j < Const, log m 
m x=- 


m 


Since 0 < R(k) <1 we have 


|r(x)| < 2 * l l R < k >l < 
m | k| <2m 


furthermore, summation by parts yields 


| sin(-|)‘R ffi (x) | - | 


y R(k)*[sln(k + i)x - sln(k - 4)x] | < 
|k|<2m 2 2 


<4+ T |R(k+l) - R(k) | • | sin(k + i)x| , 

l<|kj<2m-l 


and since R(k) are assumed to decrease monotonically 


|R ( x) I < < Const. 1 

m "|sin(i)|" 


which completes the proof 
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B. THE DECAY RATE OF THE FOURIER COEFFICIENTS-RE VISITED 


In Section 3, we concluded that the quantities 

E^(t) = H (I-S^)u^(* , t ) B 2 , satisfy for k > m, 

L (x) 

(3.14) 


the recursive inequality 


(b.l) 


V N 

E 2k^ C ^ — H — 2 maX + e 

2e/k»k ty y. 

2<T<t 


" 2ek t . E (-) 

2kV' 


In this section we complete the details for the solution of these recurrence 
relations, and obtain that for k > 2 s »m we have 


(b.2) 


E 2k (t) £ (- 


8 S E„*N 2 s+1 8 s E n »N 2 s .-s+1 ,2 t 

0 ^ t, , fi , 0 *ek t r 

-T-) * E 0 + t 1 + ~^T7 ) e - E 


e/Tc» k 


e/Tc»k 


O’ 


i.e., (3.15) holds. For s = 0, (b.2) is reduced to (3.13); now assume that 
(b.2) is valid for any k > 2 s »m: in particular, for k > 2 s+ ^*m we can 
use (b.2) with k replaced by ^ ’ Y ^ 2S * m » anc * obtain that 


8 S+1 E,,.N 2 s+1 8 S+1 E n N 2 s . — S , 2_ 

(b.3) max E (x) < ( ^° - - . -) *E n + (l + ^-) ‘e -4 ek fc .E 

t, . e/k»k e/k*k 

2<T£t 


0* 


Furthermore, we have 


(b.4) 


8 S E N 2 s+1 
E 9t- (7 ^ ^ “ — t) 


e/k*k 


8 s E n «N^ s 0 ,-s ,2 

E o + + ^ •"• 4 * *V 


Using (b.3) and (b.4) to upper bound the right hand side of (b.l) we find 
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2 „s+l„ .,2 . 

E «N 8 E »N s+1 

E„.(t> < -2 -t • ( -*n + 


2k v - o_ /tt. i.2 


2e/k*k i e/k»k 

2 


E_N“ 8 S+1 E N 2 s _.-s .2 

0 r, . 0 'i 4 ek - c .e 


2e/k * k 


(i + 


• e 


e/Ic*k 


0 


«. - 2 E k 2 t »V\* 1 . B 

+ e W 5 ° 


+ e 


-2 E k 2 t 8 V\ S - 2 .4- S sk 2 t 


.(1 + -^-y) 

e/k«k 


• e 


V 


8 S+1 E„»N 2 s+2 


The first of the four terms on the right is less than ( — - 2 ) * E f 


e/k*k 


the sum of the remaining three terms does not exceed 


8 S+1 E .N 2 S+ 1 _ 4 -s ek 2 t 

(1 + ~ 2 — ) 

e/k»k 


•e 




and hence for k > 2 S+1 *m we have 


E 2k ^ t ' ) - ^ ^ — 1 
e/Tc*k 


8 S+1 E.»N 2 s+2 8 S+1 E *N 2 s+1 ,~s .2 

0 ) •V 0 *- 7 X) 6 - E ' 


e/tc*k 


which completes the induction proof of (b.4)< 
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